Thermal Suppression of Strong Pinning 
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We study vortex pinning in layered type-II superconductors in the presence of uncorrelated disorder 
for decoupled layers. Introducing the new concept of variable-range thermal smoothing, we describe 
the interplay between strong pinning and thermal fluctuations. We discuss the appearance and 
analyze the evolution in temperature of two distinct non-linear features in the current-voltage 
characteristics. We show how the combination of layering and electromagnetic interactions leads to 
a sharp jump in the critical current for the onset of glassy response as a function of temperature. 



PACS numbers: 74.60. Ge, 74.60. Jg 

Quenched disorder in strongly layered superconduc- 
tors, such as the Bi-based high-T c compounds or the 
organic BEDT-based materials, naturally leads to the 
phenomenon of strong vortex pinning. With the mag- 
netic field directed perpendicular to the layers the vortex 
lines divide up into loosely coupled strings of "pancake" 
vortices |Q. At low fields, the absence of interactions 
between the pancake vortices allows for their free accom- 
modation to the pinning potential, thus leading to strong 
pinning of individual pancake vortices. This is in con- 
trast to the weak collective pinning situation |^], where 
the disorder potential competes with elastic forces, either 
due to tilt or shear energies. Whereas a detailed under- 
standing of the weak-pinning phenomenology has been 
developed over the past decade ||, not much progress 
has been made regarding the strong-pinning situation. 

Recent experimental and theoretical interest concen- 
trates on the low-field properties of vortex matter in 
strongly layered materials such as Bi2Sr2CaiCu2 08 (BiS- 
CCO) with particular emphasis on the effects 

of thermal fluctuations, quenched disorder, and their mu- 
tual interplay. In this letter, we present a detailed analy- 
sis of the phenomenon of strong vortex pinning in layered 
type-II superconductors and its crossover to the weak- 
pinning situation due to thermal fluctuations. We assume 
quenched point-like disorder and concentrate on the de- 
coupled limit where electromagnetic interactions between 
pancake vortices determine the behavior of vortex mat- 
ter. We discuss the presence of two step-like features in 
the current-voltage characteristics (CVC) and determine 
the evolution in temperature of the two corresponding 
critical current densities j pc (T) and j g (T) describing the 
onset of strong pinning and of glassyness, see Fig. [I]. To 
account for thermal fluctuations smoothing the disorder 
potential in the strong-pinning regime we introduce the 
new concept of variable-range thermal smoothing. We 
show how the thermal depinning of vortices proceeds in 
a sequence of steps until the usual weak collective pin- 
ning situation is recovered at high temperatures. Fur- 
thermore, we predict that the strong dispersion in the 
electromagnetic tilt modulus leads to a pronounced jump 



in the critical current density j g with increasing temper- 
ature. The main results of our analysis are summarized 
in the pinning diagram of Fig. ^| which shows the vari- 
ous pinning regimes present in the low-field region below 
the field B\ = $o/A 2 , where A is the planar London 
penetration depth and $o denotes the flux unit. In a 
forthcoming publication [p) we will map out the entire 
pinning diagram and generalize our analysis to include 
finite Josephson coupling. 




FIG. 1. Current- voltage characteristics of layered super- 
conductors in the pancake-vortex pinning regime (T < T pc ) 
exhibiting a two-step behavior. Above the critical current 
density j pc we find the usual flux-flow (FF) regime with resis- 
tivity p FF . As the current j drops below j pc , pancake vortices 
are trapped into potential wells of depth U pc and the system 
enters a second ohmic regime with a reduced resistivity due 
to thermally activated flux flow of individual pancake vor- 
tices, Ptaff ~ Pff exp(— Upc/T). At jg < jp C the motion of 
pancake vortices is inhibited by elastic forces, resulting in a 
sharp drop of the voltage (glassy response). Inset: temper- 
ature dependence of the two step features in the CVC. The 
critical current density j pc remains constant up to T = U P c, 
where the corresponding step in the CVC disappears. The 
critical current density for glassy response j g decreases with 
temperature and, at T = T pc , smoothly goes over to the crit- 
ical current density j c determined by weak collective pinning 
theory. Note that T pc 



Upc as shown in the text. 



We begin our study with the low-field/low-temperature 
limit and consider an individual vortex line oriented per- 
pendicular to the layers. The presence of point disorder 
leads to a distortion of the vortex line with a typical rel- 
ative displacement u between neighboring pancake vor- 
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tices. The optimal pinning state is determined by the 
competition between the elastic and the pinning ener- 
gies. The deformation of the vortex line on a length 
L costs an energy £ c \{u,L) ~ £i{u, k z ~ 1/L)u 2 /L. 
For purely electromagnetic interaction (uncoupled lay- 
ers) the elasticity e; takes the strongly dispersive form 
ei(u < A, k z ) = (e /2X 2 k 2 ) ln[l + \ 2 k 2 J{l + u 2 k 2 )}, with 
the line energy e = ($ /47rA) 2 . Here, we have inter- 
polated between the formulae valid for uk z > 1 M and 
uk z < 1 0]. On the other hand, adjusting to the disorder 
potential, a vortex segment of length L gains the pinning 
energy £ p ; n (u, L) ~ |£T (k) | i/c?, where d is the layer 
separation and £q(u) is the deepest minimum a pancake 
vortex can settle in within the area u 2 . This energy is 
determined by the condition u 2 J £ ° g{£)d£ ~ 1, where 
g{£) is the distribution of pinning energies, which for a 
large number of defects we assume to be Gaussian B 
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(i) 



Here, U p quantifies the disorder strength and £ is the 
planar coherence length (and also the typical distance 
between pinstates). For u 3> £ (strong pinning) each 
pancake vortex can explore many minima and one finds 



£o(u) 



In 
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< -u„ 
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Introducing the energy scale E em = e d^ 2 /A 2 , we arrive 
at the vortex free energy / per unit length, 
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Minimizing / with respect to u and L provides us with 
the optimal pinning state. For strong pinning the mini- 
mum is realized by the OD pancake-vortex configuration 
(L = d) and minimizing Eq. m|) with respect to u, we 
obtain the optimal search area fq] 
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The activation barrier for pancake- vortex motion is 

1/2 
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Comparing the Lorentz force j((fr /c)d with the pinning 
force U pc /£, we find the pancake critical current density, 



Jpc — Jp 
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Here, j p — j (U p /e d) and j a — ce /^ denotes the 
depairing current density. In the weak-pinning situation 



Up < E cm we have u g ~ £q 



-f/p, and it is ener- 



getically more favorable for the system to settle in the 
1 D pinning regime: Minimizing f(u = £, L) with respect 
to L, we find L c ~ A(A/d) 1 / 3 ( J B om /t/p) 2 / 3 > A. The 
parameter C/ p can be estimated from experiments mea- 
suring the critical current density at low B and T and is 
typically of the order of 10 K (BiSCCO). In comparison, 
the electromagnetic elastic energy E cm w 0.2 K is much 
smaller. From these estimates we conclude that we usu- 
ally encounter a strong-pinning situation with U p ^> E em 
in strongly layered high-T c material. 




FIG. 2. B-T-pinning diagram of a layered superconduc- 
tor with decoupled layers (log scales, t = T/T c , b — B/B C 2\ 
B\ = $o/A 2 ; parameters for BiSCCO: A = 1800 A, £ = 25 A, 
d = 15 A, T c = 90 K). Results for the various pinning 
regimes are shown for moderate pinning with T p — U p = 10 
K. The strong-pinning region at low temperatures T < T pc is 
divided into two parts: At T < T g thermal effects are irrele- 
vant [(0-0)D regime], while at T g < T < T pc variable-range 
thermal smoothing takes place [(0-l)D and (0-3)D regimes]. 
The dashed fat line at T — T g indicates the jump in the col- 
lective pinning length L g and the critical current density j g . 
At temperatures T > T pc and low inductions vortex segments 
of length L c > A are pinned collectively (ID regime). The 
remaining area shows the 3D (vortex-bundle) pinning regime. 

The strong pinning of individual pancake vortices into 
potential wells of depth U pc leads to a sharp drop in 
the current-voltage characteristics (CVC) at the critical 
current density j pc , see Fig. However, at finite temper- 
atures the CVC does not drop to zero, as the individual 
pancake vortices can overcome their finite pinning bar- 
riers by thermal activation, and we arrive at a second 
ohmic regime at small current densities j < j pc with a 
reduced resistivity described by thermally activated flux 
flow p TAFF - p FF exp(-C/ pc /T) [here, p FF = {B/H C2 )p n is 
the usual flux-flow resistivity, with H C2 the upper critical 
field and p n the normal-state resistivity] . Glassy response 
appears only at low current densities j < j g <C j pc , when 
the free thermal hopping of pancake vortices is inhibited 
by the elastic coupling to other pancake vortices. To de- 
termine j g we have to consider the hopping process of in- 
dividual pancake vortices. Following the usual variable- 
range-hopping (VRH) argument |l0| ], a pancake vortex 
can move freely as long as the current compensates for 
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the energy 5£(u) required to hop on to the next favorable 
state. The "minigap" 5£{u) is obtained from comparing 
neighboring favorable states: With u 2 J 1 d£g{£) = 2 we 
find 5£{u) =£ 1 -£ ~ U p [ln{u 2 f^ 2 )}- 1 / 2 . The onset of 
glassy response is determined by the condition j g {$o/c) 
u n d : 



igu, — S£(u g ), resulting in the critical current density 



3g — Jp 



1/2 



In' U > 



and the corresponding pinning energy 



Ug = 5£{Ug) ~ Up 
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As the external current j decreases below j g , single pan- 
cake vortices cannot find an appropriate final state any 
more and vortex motion involves line segments with a 
length determined by the usual laws of creep dynamics 
§,L(j)^Lg(j a /j) 5/7 (here, L g = d). 

Going over to finite temperature, we note that the sit- 
uation remains unchanged for T < U g [(0-0)D regime]. 
As T increases beyond U g , the process of variable-range 
thermal smoothing (VRS) sets in: Thermal fluctuations 
push the vortices to probe an area (u 2 ) t h > u 2 gl but elastic 
forces prevent individual pancake vortices from hopping 
to favorable states, hence segments of length L g (T) > d 
will take over the creep process. The thermal hopping of 
these segments then leads to the smoothing of the pin- 
ning energy. Hence, the temperature T g = U g defines a 
(first) thermal depinning temperature in our problem. 

We proceed with a detailed analysis of VRS: The 
mean-squared thermal displacement of a free vor- 
tex segment of length L is given by (u 2 (L,T)) t h — 
fij£(dk z /2ir)(T/£i(k z )kz), and using the dispersive elas- 
ticity from above, we find 



T 



(u 2 (L,T)) 



E, 



E KTn X 2 



d<L< X 2 /d, 
X 2 /d<L, 



T dL ^ 



[for d < L < u we should account for a log-correction 
ln[(A/£) 2 i? cm /T] guaranteeing the smooth crossover to 
u 2 as (L, T) \ (d, T g )]. The thermally smoothed pinning 
energy for a vortex segment of length L is given by 



£pm{L, T) 



; ((Ag) 2 (T)) 
N(T) 



T>Tg. 
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Here, ((A£) 2 ) = (£ 2 ) — (£) 2 denotes the fluctuation in 
the pinning energy [use Eq. ([!])], which is equal to U 2 
at T = T g and tends to U 2 /2 as T -> oo. N(T) is the 
number of available states given the search area (u 2 ) t h, 
N(T) « (u 2 )th J £ ° d£g(£), with lower and upper lim- 
its N{T g ) = 1 and N(T -> oo) = (u 2 )th/£ 2 - It is the 



suppression in the number of available states, N(T) <C 
(w 2 )th/£ 2 , which distinguishes the new VRS from the 
conventional smoothing occuring within weak collective 
pinning theory JO]. This reduced smoothing is a conse- 
quence of strong pinning and is realized within the regime 
T < T pc , where the (second) depinning temperature T pc 
is determined by the condition N(T^ 



the distribution 
entire set N(T) 



1) we find T x 



pc 



pc- 



i {u 2 ) th /2C 2 . For 
At T > T pc the 
)th/£ 2 of pinning states takes part 
in the smoothing and Eq. ( |l0| ) reproduces the standard 
expression Q, £ pin (L,T) ~ Up^/^u^^L/d. 

Within the VRS regime at low temperatures T g < T < 
T pc we use suitable numerical interpolations for the error 



function [from Jg 1 



£o+T 



d£g(£)£ n , (n = 0,2)], 
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such that £pi n (d, T) goes over to U g and U p ^£ 2 / (u 2 ) t h 
as T \ T g resp. T /" T pc . The vortex free energy / per 
unit length takes the form 
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Again, the minimum of / with respect to L defines the 
pinning length L g . Inspection of Eq. ( |l2|) reveals that 
the minimum at L = d vanishes at T = T g , implying that 
L g (T) is determined by the minimum at large lengths, 



L g {T) ~ A 



A\ 1/3 ( t^ 1/2 



exp 



T 

3T n 



< X 2 /d. (13) 



The jump from L g = d to L g (T) > X at T — T g implies 
a concurring jump in the pinning energy and the critical 
current density, 



U g (T) ~ T 



x \ 2/3 /rr \ 1/2 

A \ I lg 



TJ CXP r^)' (14) 
i(T)^»(Q 3/3 (a)' /2 »p(-^). (15) 

The jump at T g is a consequence of the strong dispersion 
in the electromagnetic line tension and persists into the 
high-field regime H 7 where it matches up with the jump 
found by Kes and Koshelev [12| in their T = analysis of 
this problem. A sharp increase (jump) in the activation 
energy U by a factor of 10 with increasing temperature 
(T > 15 K) has been found in several relaxation experi- 
ments in BiSCCO material p3[ . 

At temperatures above T pc we can ignore the under- 
lying strong pinning. Repeating the above minimiza- 
tion procedure with the usual smoothed pinning potential 
pT[, we obtain the weak-pinning results (ID regime) 
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L C (T) ~ A 
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(17) 
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where we have introduced the single-vortex depinning 
temperature . Quick inspection shows that L C (T), 
U c , and j c (T) match up with L g {T), U g (T), and the on- 
set of glassy response at j g (T) when T \ T pc . 

When the collective pinning length L C (T) exceeds 
X 2 /d the mean thermal displacement grows with in- 
creasing length and pinning becomes marginal The 
determination of L C (T) involves the calculation of the 
disorder-induced fluctuations ((u 2 (L,T))), where ((•)) 
denotes averaging over thermal fluctuations and dis- 
order. Repeating this calculation for the case of 
electromagnetic coupling, we find ((u 2 (L, T)))/£ 2 ~ 
{U p /T) 2 {L/d)\n{Ld/X 2 ). As the disorder-induced fluc- 
tuations increase beyond the thermal ones, the system 
crosses over to the pinning-dominated regime. The con- 
dition ((u\{L,T))) ~ (u 2 (L,T)) th determines the length 



and the critical current density takes the form 
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The above results properly match the previous ones at 
the crossover temperature T = T£ p . The complete tem- 
perature dependencies of the (critical) current densities 
jpa jc, an d j g are illustrated in the inset of Fig. |l|. 

As we increase the magnetic field, interactions with 
other vortices start to interfere with the accommodation 
of the individual vortices to the disorder potential. The 
regime of validity of the above results is found by compar- 
ing the tilt elastic energy E em {u/£) 2 with that of shear, 
c 66 u 2 d, with the s hear modulus given by c 66 ~ s /a 2 
and {e /X 2 )^J X/a exp(— a /X) above resp. below B\. At 
low temperatures T < T g we find that pancake vortices 
start to interact within the planes as the magnetic in- 
duction increases beyond B\, see Fig. [| At T > T g 
we have to account for the finite lengths L g (T) and 
L C (T) of the segments involved. The crossover condition 
L ~ Aexp(a /2A) produces the result 



B^-BW^ + ln- 



To 



T g <T < r pc , 
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A T 3 



dE om U 2 



T pc <T<T! p , (22) 



B 13 (T) ~ B, 
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dp 
T 



n P < t. 
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Above B 13 (T) pinning involves relaxation of vortex bun- 
dles [(0-3)D and 3D regimes]. These results apply to the 
vortex solid regime below the melting line B m (T) JL4||: 
Upon melting, both the shear and tilt moduli vanish in 
the resulting pancake-vortex gas phase, cutting off the 
ID and 3D pinning regimes at B m (T). 

Our analysis (and its generalization to higher fields 
) sheds light on two recent experiments in layered BiS- 
CCO material. The Maley analysis of the creep bar- 
rier as measured by van der Beek et al. g is consis- 
tent with a diverging barrier U(j — > 0) — > oo, rather 
than the constant barrier expected for strong pancake 
pinning in the low-T/low-B domain. The above results 
explain how individual pancake vortices couple into vor- 
tex lines exhibiting diverging barriers as the driving force 
vanishes. Second, recent local Hall-probe measurements 
of the current flow in BiSCCO crystals show a sharp, 
roughly field-independent onset of bulk pinning as the 
temperature decreases below T ss 40 K M. This exper- 
imental finding is in agreement with the appearance of 
strong non-linearities in the CVC below the temperature 
T pc , see Fig. [j]. The definite identification of this pinning 
onset with either j pc or j g requires a detailed amplitude 
and frequency analysis of the experimental feature. 
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